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FOREWORD 


This  report  was  prepared  "by  Midwest  Research  Institute  under  USAF 
Contract  No.  AF  33{616)-7058.  The  contract  was  initiated  under  Project  No. 
7360,  "The  Chemistry  and  Physics  of  Materials,"  Task  No.  736002,  "Nondestructive 
Methods."  The  work  was  administered  under  the  Directorate  of  Materials  and 
Processes,  Deputy  for  Technology,  Aeronautical  Systems  Division,  with  Mr-  R.  R. 
Rovand  acting  as  project  engineer. 

This  report  covers  work  conducted  from  1  February  1961  to  31  January 

1962. 


The  work  was  performed  under  the  direction  of  Mr.  Fred  Hollins,  Jr. 
Personnel  directly  involved  in  prosecution  of  the  research  have  been  Messrs. 
Jerry  Jones,  Donald  Kobett,  Neil  Abbott,  James  Gravitt,  Paul  Gutshall,  and 
Phillip  James. 


ABSTRACT 


The  study  of  stress- induced  birefringence  has  been  continued  in  both 
polycrystalline  and  single  crystal  experiments.  The  effect  is  explained  on  the 
basis  of  nonlinear  elasticity  theory.  Experiments  indicate  that  dislocation 
activity  does  not  strongly  affect  results  in  polyciystalline  specimens,  however, 
a  pronounced  influence  may  be  observed  in  single  crystals. 

Nonlinear  elasticity  theory  has  been  used  to  investigate  the  inter¬ 
action  of  two  intersecting,  plane,  elastic  waves  in  a  homogeneous,  isotropic 
medium.  A  criterion  for  the  occurirence  of  a  strong  scattered  wave  has  been 
derived.  The  criterion  is  formulated  as  a  relationship  between  the  second 
order  elastic  constants  of  the  material,  the  angle  between  the  intersecting 
wave  vectors,  and  the  ratio  of  primary  wave  frequencies.  The  amplitude  of  the 
scattered  wave  is  found  to  be  proportional  to  the  volume  of  interaction  and 
dependent  on  the  third  order  elastic  constants  of  the  material.  Preliminary 
efforts  to  experimentally  verify  the  theoretical  predictions  are  described. 


This  report  has  been  reviewed  euid  is  approved. 


W.  J.  TRAPP  ' 

Chief,  Strength  and  Cynamics  Branch 
'fetals  and  Ceramics  Laboratory 
Directorate  of  Materials  and  Processes 
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I .  INTRODUCTION 


In  recent  years  considerable  effort  nas  been  directed  toward  the  de¬ 
tection  and  measurement  of  residueil  stre^-es.  Most  of  the  methods  that  have 
been  successfully  developed  are  destructive  or  semidestructive  in  nature.  These 
methods,  in  general,  are  based  on  the  fact  that  relieving  the  stresses  will  al¬ 
ter  the  dimensions  of  the  body.  Successive  sections  of  the  test  piece  are  re¬ 
moved  by  chemical  or  mechanicad  means  and  the  strain  relaxation  that  occurs  in 
the  remaining  material  is  measured.  The  stress  distribution  is  then  calculated 
from  the  strain  distribution  found  in  the  metal  specimen. 

The  only  nondestructive  method  that  is  successfully  used  in  measur¬ 
ing  residual  stresses  is  one  which  utilizes  X-ray  diffraction  techniques.  Even 
this  method  has  a  number  of  disadvantages.  Good  accuracy  is  obtained  only  with 
specimens  that  yield  sharp  diffraction  lines.  Quenched  steels  or  cold  worked 
metals  give  diffuse  lines  that  produce  very  large  errors .  An  even  greater  lim¬ 
itation  is  that  only  surface  stresses  can  be  investigated  using  this  method. 

During  the  past  two  years  Midwest  Research  Institute  has  been  engaged 
in  a  study  of  the  stress -dependent  aspects  of  ultrasonic  propagation  in  solid 
materials .  It  is  hoped  that  the  results  of  this  study  may  facilitate  the  ap¬ 
plication  of  ultrasonics  to  specific  residual  stress  problems .  An  earlier  re- 
porti/  covers  the  first  year's  work  and  snouid  be  consulted  for  experimented,  de¬ 
tails  that  have  been  intentionally  oniitted  in  the  following  discussion. 


II.  STRESS-INDUCED  BIREFRINGENCE 


A.  General  Discussion 


It  has  been  recognized  for  some  time  that  the  velocity  and  attenua¬ 
tion  of  ultrasonic  waves  traveling  throu^  solid  materled.s  are  usually  stress- 
dependent.  The  exact  nature  of  these  stress-dependent  changes  is  not  completely 
understood  but  recent  advances  in  dislocation  theory  and  finite  elasticity  can 
expleiin  many  of  the  experimental  observations . 

A  property  that  we  have  studied  in  several  solids  is  the  stress-in¬ 
duced  double  refraction  of  megacycle  shear  waves.  In  a  rectangular  coordinate 
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system^  let  us  consider  a  shear  wave  traveling  in  the  x-directlon.  If  the  ma- 
terisd  is  isotropic,  we  will  find  that  the  shear  wave  velocity  is  independent 
of  the  direction  of  particle  motion.  However,  if  we  remove  the  isotropy  by 
applying  a  tensile  stress  in  the  y-direction,  the  velocity  is  found  to  vary 
with  the  direction  of  particle  motion.  Under  these  conditions  it  is  impossible 
to  propagate  a  plane  polarized  shear  wave  except  when  the  particle  motion  is 
either  parallel  or  perpendicular  to  the  applied  stress.  Referring  to  the  above 
two  velocities  as  Vgx  and  Vs2  respectively,  we  can  define  a  fractional  ve¬ 
locity  difference,  (dv/v)  =  (V3i-Vs2)/V  ,  where  V  is  an  average  of  Vgl  and 
Vgg.  The  value  of  Av/v  can  be  determined  experimentally  using  a  simple  pulse 
echo  technique  .i'  It  ha.s  been  well  established  that  the  Av/v  values  are 
stress-dependent  and,  in  fact,  vary  linearly  with  the  stress  level.  The  cor¬ 
relation  between  theory  and  observation  is  described  in  the  following  two  sec¬ 
tions  . 


B.  Predictions  Based  on  Theory  of  Finite  Strain 

The  velocities  of  longitudinal  and  shear  waves  in  an  isotropic,  homo 
geneous  medium  are  usually  given  by  the  following  expressions: 


where  X  and  u  are  the  second  order  Lame  constants  and  p  is  the  density. 
These  expressions  are  derived  on  the  basis  of  infinitesimal  deformations.  They 
do  not  predict  a  stress -dependent  velocity  so  long  as  the  density  remains  con¬ 
stant.  However,  by  using  the  finite  strain  method  of  Mumaghan,  2/  Hughes  and 
Kelly^/  calculated  the  velocities  of  ultrasonic  waves  in  isotropic  media  and 
obtained  the  stress -dependent  expression  shown  below: 

(a)  For  hydrostatic  pressurej  p 


Longitudinal  wave  Vl  =  (2ui+^)  “  ^  [lOp+7X+6£+4i0 


(2) 


Shear  wave 


PoVi  =  U  -  [3(2u+x)+3m  -  |] 


(3) 
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(b) 


(c) 


Uniaxial  Pressure, 


Longitudinal  wave 


Slieur  wave 


Uniaxial  Pressure, 
direction 


P  ,  and  wave  propagation  both  in  x-direction 
PqVl  =  (2ii+X)  -  ^  (l0a+4X+4m)+X+2jej 

i’o'i  “  I*  - 1(  ^ 

P  ,  in  x-direction  \cith  wave  propagation  in  y- 


(4) 

(5) 


Longitudinal  wave 


Shear  wave  with 
peurticle  motion 
in  x-direction 


Po'^L  =  (2u+X)  -  L-  [2/  -  ^  (24+X-Hn)] 

Ooil  =  “  -  ^  ”] 


(6) 


(7) 


Shear  wave  with 

particle  motion  PoVio  =  n  -  rm-2X  -  ^  nl  (8) 

in  2 -direction  ^  ^  J 


In  these  expressions,  Pq  is  the  density  of  the  undefomed  medium,  K  is  the 
bulk  modulus  and  the  terms  i  ,  m  ,  and  n  are  referred  to  as  third  order 
elastic  constants .  The  latter  give  an  indication  of  the  nonlinearity  in  the 
elastic  behavior  of  a  material.  It  is  easily  seen  that  when  the  hydrostatic 
pressure  or  uniaxial  stress  is  zero  the  expressions  all  reduce  to  the  two  ex¬ 
pressions  given  in  (1) . 

The  three  third  order  constants  that  appear  in  expressions  (2)-(0) 
can,  of  course,  be  determined  experimentally  by  accurately  measuring  the  ve¬ 
locity  changes  produced  in  at  least  three  independent  cases.  It  is  possible 
to  independently  determine  the  third  order  constant,  n  ,  by  a  simple  technique 
that  does  not  necessitate  the  accurate  measurement  of  absolute  velocity.  When 
shear  waves  are  sent  through  a  specimen  in  a  direction  that  is  perpendicular 
to  the  direction  of  an  applied  uniaxial  stress,  expressions  (7)  and  (8)  in¬ 
dicate  that  the  velocity  depends  upon  the  particle  motion  direction.  The  in¬ 
itially  isotropic  material  becones  birefringent .  Pure  shear  waves  can  be 
propagated  in  this  case  only  when  the  particle  motion  is  either  parallel  or 
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perpendicular  to  the  applied  Aoaa.  If  we  subtract  expression  (0)  from  (7)  and 
remember  that  K  =  l/3(2^l+3X)  ,  we  get 


OoVfi  "  [4p+n  j  . 


(9) 


Since  +  V^2  —  ^^,3  ,  we  then  have 


^32  - 


ePo^oli 


[4n+n] 


(10) 


or  expressed  as  a  fractional  difference 


’^S1"''^S2 


=  dV/V  ^ 


|_4M,+nJ 


^  [4^+n]  (11) 

0)J, 


From  (11)  we  see  that  the  fractional  difference  between  and  Vg2  is  a 

linear  function  of  the  uniaxial  stress^  P  .  Note  also  that  experimental  deter¬ 
mination  of  iv/v  at  a  given  stress  is  sufficient  to  evaluate  the  third  order 
constant,  n  . 

The  velocity  of  ultr£-sonic  waves  in  single  crystals  is  also  stress- 
dependent,  but  the  velocity  expressions  are  somewhat  more  involved  due  to  the 
loss  in  symmetry.  Seeger  and  Buckl/  have  derived  the  differential  equations 
for  wave  propagation  in  cubic  crystals  that  have  superimposed  finite  stresses. 
There  are  six  third  order  elastic  constants  for  such  crystals  and  Bateman 
et  al.,^/  have  recently  used  ultrasonic  techniques  to  evaluate  all  of  these 
constants  for  germanium.  The  velocity  expression  given  in  Ref.  5  for  a  shear 
wave  traveling  in  the  joo^  direction  with  particle  motion  parallel  to  a  [lid 
pressure  is  as  follows: 

0V||  =  Cl+7+a-e)C44  +[(tt+B)/4]  0^44  +  [(a+e+2Y)/4]  Cj^gg 

+  [(a-P)/-^  C456  .  (12) 
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For  a  shear  wave  traveling  in  the  same  direction  but  with  particle  motion  per¬ 
pendicular  to  [ll^  ,  the  expression  is 


pVj^  =  (l+Y+!?-a)C44  +  [(a+0)/4]  C144  +  [(a-*-p+2Y)/4] 

+  [(p-a)/4]  C456  .  (13) 


In  the  above  expressions  a  ,  P  ,  and  Y  are  all  linear  functions  of  the  ap¬ 
plied  stress  and  the  C4 terms  are  the  third-order  moduli •  By  using  the 
approximation  (V„  +  Vj_  )  2:  2V  ,  the  fractional  velocity  difference  for  shear 

wave  propagation  in  the  [iCX^  direction  becomes 


^^^44'^C456^  • 

^44 


(14) 


Here  again  we  see  that  the  fractional  velocity  difference  is  a  linear  function 
of  the  uniaxifij.  stress  P  . 

It  should  be  emphasized  at  this  point  that  the  above  expressions^ 
(2)-(8)  and  (12)-(13),  were  derived  for  a  perfectly  elastic  material.  For  most 
real  macerials,  an  applied  stress  will  produce  some  anelastic  strain  in  addi¬ 
tion  to  the  elastic  strsdn .  Granato  and  Lucke^zZ/  have  developed  a  theory  of 
energy  losses  and  modulus  changes  that  are  due  to  dislocation  damping.  The 
theory  is  based  on  the  model  of  a  dislocation  loop  oscillating  under  the  in¬ 
fluence  of  an  applied  stress.  Without  going  into  the  details  of  this  theory, 
it  is  sufficient  to  state  that  the  fractional  velocity  change  due  to  the  dis¬ 
locations  is  proportional  to  the  dislocation  density  and  the  second  power  of 
the  dislocation  loop  length.  There  is  considerable  experimental  evidence  that 
supports  the  Granato-Lucke  theory  and  there  is  no  question  that  dislocations 
do  affect  the  propagation  of  ultrasonic  waves.  Additional  discussion  of  this 
subject  may  be  found  in  the  following  section. 


C .  Experimental  Results 

i,ne  puxse-echo  techniquei/  has  been  used  to  measure  stress-induced 
•hanges  of  (AV/v)  in  a  variety  of  materials.  Figure  1  shows  some  typical 
urves  <’nd  confirms  the  linear  relationship  between  (AV/v)  and  stress  level. 
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Figure  1  -  Stress- Induced  Birefringence  for  Several  Materials 


It  was  mentioned  previously  that  the  (dv/v)  measurements  might  be  used  to 
evaluate  the  third  order  elastic  constant,  n  .  Vfe  see  from  Eq.  (ll)  that  a  plot 
of  (fiv/v)^  versus  stress  should  yield  a  straight  line  with  a  qlope  of 
“(4|j,+n)/8y,  .  The  slopes  of  the  curves  shown  in  Figure  1  have,  therefore, 

been  used  to  calculate  n  for  each  material  and  the  results  are  summarized 
in  Table  I. 


TABLE  I 


Material  n  (psl) 


6061-T6  Aj-umlnum 

-(45.2  t  1.8) 

X  10® 

UOO-F  Aluminum 

“(69.2  -  5.5) 

X  106 

C1018  Stet^ 

-(83.1  t  4.4) 

X  106 

Fused  Silica 

-(33.2  1  0.9l 

X  106 

It  would  be  interesting  to  conq)are  the  values  in  Table  I  with  those  obtained 
by  other  workers  but  we  have  been  unable  to  find  any  third  order  data  on  these 
specific  materials. 

It  has  been  previously  reported^/  that  (AV/v)  for  C1018  steel 
varies  linearly  with  elastic  strain  but  is  relatively  insensitive  to  plastic 
strain.  This  result  suggests  that  the  dislocation  activlty,which  occurs  in 
the  plastic  range,  must  affect  Vg3_  and  Vg2  nearly  equal  amounts  *  Thus, 
(Vsi-Vs2)  remains  essentially  unchanged.  The  variation  of  (AV/v)  during 
elastic  and  plastic  defomation  has  more  recently  been  studied  in  specimens 
of  1100-F  aluminimi.  The  results  have  been  very  similar  to  those  for  steel. 
Figure  2  presents  data  for  a  specimen  of  UOO-F  eiluminum  both  before  and  after 
yielding.  The  first  data  were  taken  during  the  initial  loading  cycle.  The 
specimen  was  then  loaded  beyond  the  yield  point  until  a  strain  of  approximately 
per  cent  was  reached  and  then  the  second  set  of  data  was  taken.  The  maxi¬ 
mum-minimum  type  pulse  echo  pautems,  which  are  used  in  evaluating  (AV/v)  , 
change  continuously  as  the  load  varies  but  they  do  not  appear  to  vary  with  time 
at  constant  load.  However,  if  the  transducer  is  rotated  so  that  the  particle 
motion  is  parallel  or  perpendicular  to  the  cctnpression,  an  exponential  decay 
pattern  is  obtained  and  the  attenuation  is  very  time -dependent  innnediately  fol¬ 
lowing  load  changes.  The  time -dependent  variation  in  attenuation  is  apparently 
due  to  dislocations  that  break  away  from  impurities  and/or  vacancies  when  the 
stress  is  varied.  These  point  defects  then  diffuse  through  the  lattice  and 
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Figure  2  -  Stress-Induced  Birefringence  in 
1100-F  Aluminum 
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repin  the  dislocations  as  the  load  is  held  constant.  In  contrast,  the  lack  of 
time  dependence  in  the  stress-induced  birefringence  together  with  the  agreement 
between  pre-yield  and  post-yield  values  provide  a  strong  indication  that  dis¬ 
locations  have  little  effect  on  birefringence  measurements  in  polycrystalline 
steel  or  aluminum* 

Me  have  also  studied  stress-induced  birefringence  in  single  crystals 
of  high  purity  aluminum  and  sodium  chloride .  The  results  have  been  similar  to 
those  obtained  for  polycrystalline  specimens  but  somewhat  less  reproducible. 
There  is  an  additional  effect  which  one  may  encounter  in  working  with  single 
crystals  and  that  is  internal  conical  refraction.  This  effect  occurs  when  shear 
waves  travel  through  a  crystal  in  a  direction  that  is  not  exactly  parallel  to 
a  specific  crystal  direction.  Pure  shear  waves,  for  example,  can  travel  through 
a  cubic  crystal  in  the  [lOO]  direction  and  the  velocity  is  independent  of  the 
polarization  angle.  If,  howler,  the  propagation  direction  is  slightly  mis- 
oriented  with  respect  to  [lO^  ,  then  internal  conical  refraction  may  occur. 
V/atennan  and  TeutonicoS/  have  shown  that  this  effect  can  produce  pulse-echo 
patterns  exactly  Uke  those  used  to  evaluate  stress -induced  anistropy.  Precise 
orientation  of  the  single  crystal  is  very  important  then  in  shear  wave  studies, 
particularly  at  high  frequencies.  The  technique  of  Walker  and  co-workers^/ 
was  used  to  orient  and  grind  the  mechanical  faces  of  the  specimen  parallel  to 
a  desired  crystal  face.  Subsequent  hand  lapping  was  necessary  to  attedn  the 
even  more  stringent  requirements  on  parallelism  between  opposite  surfaces .  For 
pulse.echo  work  at  frequencies  between  10  and  100  me .  the  parallelism  should 
>'e  maintained  to  within  0.00005  in/in. 

In  addition  to  the  problems  outlined  above,  the  single  crystal  experi¬ 
ments  were  influenced  strongly  by  dislocation  activity.  In  both  aluminum  and 
sodium  chloride  the  shear  waves  were  sent  throu^  the  crystal  along  (l0(^ 
while  compressive  loads  were  applied  in  the  [ooj  direction.  As  the  external 
stress  was  applied  to  aluminum  the  (Av/v)  values  were  observed  to  be  both 
stress  dependent  and  time  depen«^ent.  Upon  loading  the  specimen  (Av/v)  first 
Increased  and  then  part  of  the  increase  was  lost  as  time  progressed  at  constant 
load.  The  equilibrium  value,  reached  after  a  few  minutes,  was  proportional  to 
the  applied  load.  V/hen  the  load  was  removed,  (Av/v)  mcxnentarily  increased 
again,  but  then  decreased  with  time  to  the  no  load  value.  The  equilibrium 
value  after  loading,  however,  was  usually  higher  (indicating  increased  aniso¬ 
tropy)  than  nrlor  to  loading. 

In  similar  experiments  on  sodium  chloride  crystals  the  (AV/v)  valuer 
.rere  more  reproducible  and  less  time  dependent  than  for  aluminum.  The  varia¬ 
tion  is  probably  due  to  the  different  slip  systems  found  in  the  two  materials . 
rhe  primary  slip  system  for  aluminum  is  along  (ill)  planes  whereas  sodium 
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chloride  slips  along  (llO)  planes .  For  pure  shear  waves  traveling  in  a 
ClO^  direction  there  is  no  resolved  shear  stress  across  (110)  planes  but 
there  is  a  resolved  stress  across  (ill)  planes.  Thus,  we  see  that  a  [lO^ 
shear  wave  could  interact  with  the  many  (Ul)  dislocations  in  aluminum  and 
yet  exhibit  no  Interaction  with  the  (110)  dislocations  in  sodium  chloride . 

A  change  in  (iv/v)  due  to  dislocations  ’..•ould  only  arise,  however,  from  an 
anisotropic  change  in  the  distribution  of  dislocations.  Such  anisotropy  is 
highly  probable  in  single  crystals  where  slight  changes  in  the  direction  of 
stress  application  can  activate  more  dislocation  activity  on  some  primary  slip 
systems  than  on  others . 


Ill .  ULTRASOMIC  BEAM  INTERACTION  IN  SOLIDS 


A .  Theory 


The  nondestructive  measurement  and  analysis  of  a  three-dimensional 
stress  distribution  is  a  difficult  problem  to  contemplate.  It  is  obvious  that 
to  succeed  in  this  task  one  must  monitor  some  property  that  is  stress  sensi¬ 
tive  and  must  also  localize  the  volume  elements  within  the  test  specimen  \/hich 
affect  the  property  being  measured.  Ultrasonic  \/aves  show  some  promise  of 
being  successfully  applied  to  this  problem.  In  addition  to  certain  stress- 
dependent  aspects,  ultrasonic  vraves  can  be  collimated  into  well  defined  beams. 
The  intersection  of  two  collimated  beams  can  be  used  to  define  a  reasonably^ 
small  volume  element  within  a  relatively  large  metal  specimen.  Beam  inter¬ 
action  or  scattering  at  the  point  of  interaction  might  be  related  to  stresses 
at  the  point  of  intersection,  l/hether  or  not  the  beam  interaction  is  stress 
sensitive  is  secondary,  of  course,  to  whether  there  is  a  detectable  beam  inter¬ 
action  at  all. 

If  one  considers  beam  interaction  based  on  linear  elasticity  theory 
it  is  easily  shown  that  the  separate  waves  can  each  exist  independently  of  the 
other,  i.e.,  they  do  not  interact.  However,  if  the  basic  expression  of  elastic 
energy  is  extended  to  terms  cubic  in  strains  the  equations  of  particle  motion 
then  contain  quadratic  terms  in  the  displacements.  The  independent  nature  of 
the  intersecting  waves  no  longer  holds  and  the  theory  admits  the  possibility 
of  interaction  between  waves.  An  analysis  of  wave  interaction  in  nonlinear 
solids  has  been  completed  during  the  past  year  and  the  results  are  encourag¬ 
ing.  A  rather  detailed  description  of  the  analysis  is  given  in  the  Appendix 
so  the  results  will  be  only  briefly  summarised  here. 

The  analysis  considers  the  interaction  of  two  waves  in  a  homogenous, 

'•  sotropic  solid .  The  results  indicate  that  for  a  given  material  there  does 
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exist  an  angle  of  intersection  between  two  waves  of  frequencies  ct>j^  and 
\rtiich  produces  a  "resonant"  interaction.  Under  resonant  conditions  the  analysis 
indicates  that  scattered  waves  (with  frequencies  i  (Ug)  originate  at  the 
"point"  of  intersection.  The  analysis  further  predicts  the  direction  of  par¬ 
ticle  motion  as  well  as  the  propagation  direction  of  the  scattered  wave.  The 
resonant  angles  of  intersection  and  the  directions  of  propagation  for  the 
scattered  waves  are  summarized  in  Table  II  of  Appendix  1,  p.  26. 

Some  calculations  have  also  been  made  on  the  amplitude  of  the  scat¬ 
tered  wave.  The  amplitude  is  found  to  be  proportional  to  the  volume  of  inter¬ 
action  and  also  dependent  on  the  third  order  elastic  constants  of  the  material. 

A  numerical  example,  based  on  data  for  polystyrene,  indicates  that  the  dis¬ 
placement  amplitude  of  the  scattered  wave  is  not  beyond  the  realm  of  detection. 
Nevertheless,  the  intensity  of  this  wave  will  undoubtedly  be  stnol i  and  optimum 
conditions  must  be  obtained  if  detection  is  to  be  expected. 


B.  Preliminary  Experiments 

Experimental  efforts  to  verify  the  wave  interaction  predictions  are 
still  in  the  preliminary  stages,  but  a  description  of  the  initial  experiment 
is  included  here.  The  only  case  that  has  been  examined  experimentally  is  the 
interaction  of  two  transverse  waves.  From  Table  II  of  Appendix  I,  we  see  that 
the  scattered  wave  in  this  case  has  a  longitudinal  mode  and  a  frequency  of 
+  ujg)  .  The  table  eilso  indicates  that  the  angle  for  resonant  interaction 
is 

=  arc  cos  (a^+l)/2^ 

and  thax  the  direction  of  the  scattered  wave  is  .  We  see,  then,  that 

the  resonant  angle  depends  only  on  the  transverse  to  longitudinal  velocity 
ratio,  C  ,  and  the  frequency  ratio,  a  ,  of  the  two  primary  waves.  The  di rec¬ 
kon  of  the  scattered  wave  is  determined  by  the  propagation  vectors,  and 
K2  ,  of  che  primary  waves .  It  is  further  shown,  in  the  numerical  example  of 
Appendix  I,  that  the  amplitude  of  the  scattered  wave  is  a  maximum  when  = 

.  Under  these  conditions,  a  wave  of  frequency  2a  is  predicted  to  emerge 
in  a  direction  that  bisects  the  angle  between  the  two  primary  waves.  The  anal¬ 
ysis  also  indicates  that  the  two  primary  transverse  waves  should  both  be  po¬ 
larized  either  parallel  or  perpendicular  to  the  (i^,  i^)  plane. 


1 

> 

J 


11 


A  schematic  of  the  interaction  experiment  descrihed  above  is  shown 
in  Figure  3.  The  specimen  was  machined  with  the  shape  of  an  isosceles  tri¬ 
angle  -  the  angle  between  the  equal  sides  being  cut  so  the  two  primary  waves 
would  intersect  at  the  resonance  angle,  ^  •  The  scattered  wave  should  then 
travel  noimal  to  the  third  side  of  the  triangular  specimen.  The  two  transducers 
used  to  generate  the  primary  waves  were  both  5  me.,  AC-cut  quartz  crystals. 

The  receiving  transducer  was  a  10-mc,x-cut  crystal.  The  figure  depicts  the 
ultrasonic  wave  packets  after  they  have  passed  throu^i  the  volume  of  inter¬ 
action  and  the  scattered  wave  is  Just  arriving  at  the  receiver  crystal.  The 
electrical  signal  from  the  crystal  is  fed  to  an  amplifier  tuned  for  10  me., 
and  the  amplified  signal  drives  the  vertical  deflection  plates  of  a  cathode 
ray  indicator. 

The  transmitting  crystals  were  first  excited  simultaneously  with  a 
single  5-mc.  signal  from  one  pulse  generator.  Experiments  were  also  performed 
using  two  pulse  generators.  The  two  generators  were  triggered  at  the  same  time 
but  the  frequency  of  one  coxild  be  vetried  sli^tly.  This  procedure  was  used  to 
possibly  compensate  for  slight  errors  in  both  theory  and  experimental  arrange¬ 
ment  that  could  cause  the  angle  i  to  be  in  error.  The  compensation,  of 
course,  was  expected  from  changes  in  the  frequency  ratio  («3i/a^)  which  affectc 
the  angle,  ^  . 


Judicious  placement  of  the  three  crystals  facilitates  the  detection 
of  the  scattered  wave  by  making  it  possible  to  distinguish  the  desired  signal 
from  multiple  reflections  of  the  primary  waves .  Frexn  simple  geometrical  con¬ 
siderations,  it  is  easy  to  calculate  the  time  of  travel  necessary  to  bring  the 
primary  transverse  waves  to  the  intersection  zone  and  the  additional  time  re¬ 
quired  for  the  scattered  longitudinal  wave  to  travel  to  the  receiver  crystal. 
It  is  then  only  necessary  to  use  a  scope  with  a  calibrated  sweep  to  localize 
the  desired  signal.  Multiple  reflections  of  primary  waves  will  arrive  at  the 
receiver  crystal  much  later  than  the  scattered  wave . 

The  interaction  of  two  transverse  waves  has  thus  far  been  studied 
only  in  fused  silica  and  polycrystalline  magnesium.  We  have  not  yet  been  able 
to  detect  the  predicted  longitudinal  waves,  but  several  changes  are  being  made 
to  optimize  variables  and  increase  sensitivity.  An  increase  in  the  amplitude 
of  the  primary  waves  and  additional  amplification  of  the  received  signals  are 
possible.  Different  specimen  materials  will  be  investigated  to  optimize  non¬ 
linear  parameters.  It  is  also  desirable  to  operate  at  higher  frequencies  be¬ 
cause  the  analysis  indicates  the  amplitude  of  the  scattered  wave  is  dependent 
upon  the  third  power  of  the  primary  wave  frequency.  This  factor,  of  course, 
must  be  weighed  against  the  increase  in  attenuation  that  occurs  at  the  higher 
frequencies . 
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Figure  3  -  Schematic  of  Beam  Interaction  Experiment 


Apparatus  is  also  Ijeing  assembled  for  the  study  oi  Interactiuu  ue- 
tween  one  transverse  wave  and  one  longitudinal  wave.  This  case  has  one  advan* 
tage  over  the  transverse -transverse  case  in  that  the  longitudinal  wave  can  be 
coupled  into  solid  specimen  through  a  liquid.  This  makes  it  much  easier  to 
continuously  vary  the  interaction  angle,  ^  ,  and  thus  optimize  the  resonant 
condition.  A  complete  investigation  of  the  interaction  phenomena  will  include 
numerous  experimental  arrangements  using  both  pulsed  and  continuous  wave  tech¬ 
niques  . 


IV.  sumksci 


The  study  of  stress -induced  birefringence  has  been  continued  in  both 
polycrystedline  and  single  crystal  experiments .  The  effect  is  explained  on  the 
basis  of  nonlinear  elasticity  theory.  Experiments  indicate  that  dislocation 
activity  does  not  strongly  affect  results  in  poly crystalline  specimens,  however, 
•\  pronounced  Influence  may  be  observed  in  single  crystals. 

Nonlinear  elasticity  theory  has  also  been  used  to  investigate  the 
interaction  of  two  intersecting,  plane,  elastic  waves  in  a  homogeneous,  iso¬ 
tropic  medium.  A  criterion  for  the  occurrence  of  a  strong  scattered  wave  has 
been  derived.  The  criterion  is  foimulated  as  a  relationship  between  the  first- 
order  elastic  constants  of  the  material,  the  angle  between  the  intersecting 
wave  vectors,  eind  the  ratio  of  primal^'  wave  frequencies.  The  amplitude  of  the 
scattered  wave  is  found  to  be  proportional  to  the  volume  of  interaction  and 
dependent  on  the  third  order  elastic  constants  of  the  material.  Preliminary 
efforts  to  experimentally  verify  the  theoretical  predictions  are  described. 
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APPENDIX  I 


THE  INTERACTION  OF  ELASTIC  WAVES  IN  AN  ISOTROPIC  SQT.TD 

I.  INTRODUCTION 


In  the  linear  theory  of  elasticity  two  elastic  waves  do  not  inxeract. 
The  equations  of  motion  are  linear  and  therefore  the  principle  of  superposition 
•  Any  solution  of  the  equations  of  motion  can  be  written  as  a  linear  com¬ 
bination  of  monochromatic  waves.  The  linear  theory  of  elasticity  results  from 
assuming  the  elastic  energy  to  be  quadratic  in  the  particle  displacements.  If 
terms  cubic  in  the  particle  displacements  are  included  in  the  elastic  energy, 
the  equations  of  motion  become  nonlinear- J-O/  This  nonlineeurity  gives  rise  to 
sn  effective  interaction  between  two  plane  elastic  waves  ■vrtiich  can  produce 
scattering.  The  scattering  of  two  collimated,  monochromatic,  plane  waves  in 
infinite  isotropic  solid  is  considered  in  the  present  paper. 


II.  THEORY 


When  terms  cubic  in  the  particle  displacements  are  included  in  the 
elastic  energy,  the  resultant  nonllneju:  equatlcxis  of  motion  for  an  isotropic 
solid  arelO/ 


2  P 

a  vb 

ax,8x, 


t  . 


V^dxjjdxjj  dxj  dXj^dX]^  dXjj  dxj^  y 

+  (K+h/S+A/a+B)/^  — 

\dxjdxj^  dxjj  / 


VdXfcdXj^  dx^y  dXi  dx^dxj^  dx^y 


'  \dxidxk  dx^y 


(I-l) 


IS 


where  =  density  of  the  xande formed  solid 

component  of  the  particle  displacement 
K  =  compression  modulus 
^  s  shear  modulus. 

A  ,  B  ,  and  C  are  third  order  elastic  constants,  i.e.,  they  are  the  coef¬ 
ficients  of  cubic  strain  terms  in  the  elastic  energy.  Subscripts  appearing 
twice  In  a  single  tern  indicate  summation  over  the  values  1,  2,  5,  The  terms 
xi  ,  X2  ,  and  X3  are  rectangular  coordinates. 

The  left  side  of  Eq«  (l-l)  is  lineaur  in  iS  while  the  right  side  is 
quadratic.  In  practical  applications  the  displacement  vector  u  is  small  and 
therefore  the  right  side  of  Eq.  (l-l)  will  be  small  compared  to  tlie  left.  To 
solve  Eq.  (l-l)  we  set  =  u°+u®  where  vfi  is  the  solution  vrtien  the  right 
side  of  Eq.  (l-l)  is  zero,  and  tfs  is  a  presumably  small  correction  arising 
from  the  right-hand  side.  Up  being  a  solution  to  the  linear  equation  consists 
of  a  superposition  of  monochromatic  waves.  Since  we  are  interested  in  the 
mutual  scattering  of  two  waves  we  take 

Uq  *  Aq  ®o  ■  (l“2) 


In  Eq.  (l-2)  the  amplitude  vector  may  be  chosen  either  parallel  or  perpendicu¬ 
lar  to  the  wave  vector  It  ,  that  is,  the  primary  waves  are  either  longitudinally 
or  transversely  polarized.  For  a  transverse  wave  (a  =  C-^k  where 

^  1 
For  a  longitudinal  wave  uj  =  C^k  where  (K+4  |j,/3)/p^j^. 


Now  to  solve  Eq*  (l-l)  we  substitute  u  =  on  the  left  side 
and  'i  ^  "iP  on  the  right  side*  This  should  be  a  good  approximation  if  Ug 
is  small  compared  to  .  Since  tP  is  a  solution  of  the  linear  equation 
it  disappears  from  the  left  side  and  we  obtain 


d 


dx]^dx]^ 


(K+  ^i/3) 


dx^axi 


=  Pi 


(1-3) 
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\rtiere  the  vector  J  is  determined  by  putting  in  the  right  side  of  Eq. 

^*1*  (l”2)  for  do  ,  J  will  involve  a  sum  of  products  of 
two- monochromatic  waves.  Some  of  the  terms  will  represent  the  interaction  of 
a  primary  wave  with  itself,  ^is  interaction  has  been  treated  be fore ^O/  and 
we  shall  not  include  it  here.  Including  only  those  terms  representing  inter¬ 
actions  between  the  primary  waves  we  find  for  $  after  some  tedious  manipula¬ 
tion 


~  sin 

+  I"  sin  j(oj^-a2)t-(tj^-t2)  •  r  |  (I-4) 


where 


I*  =  -i  W4)((VBo)(k2.4)ki  i 


+  (Bo-ki)(i^-^)Ao  ±  (Ao-4)(ki.ki)Bo^(Ao.i^)(k3^.k2)Bo 
±  2(BQ*kj^)(k3^*i^)AQ  j-  ^  (K+4/3  +  A/4  +B)  {(Ao-Bo)(ki-i^  )i^ 

-  (Ao*®o)(^1*^)^1  ^  (Bo*^)(^*^)^o  -  (Ao*ki)(ki*^)Bo| 

-  i  (A/4+B)[(AQ.k2)(Bo-k2)ki  ±  (4*ki)(Bo-ki)^  XBo-k^)!^ 

*  (Ao-4)(Bo*^l)kl}-  i  (B+2C)  |(Ao.ki)(Bo-k2)k2 

+  (Ao-ki)(Vk2)ki  j 


Since  J  is  a  known  function  of  t  and  t  ,  Eq.  (1-3)  is  now  a  Unear  in¬ 
homogeneous  eqxjation  for  the  scattered  wave  tfs  ^  „ith  $  acting  as  a  source 
term.  Instead  of  writing  Eq.  (l-3)  in  c(mq>onent  form  we  prefer  to  use  vector 
notation  euad  write 
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(r,t)  -  C^7  [7.u®{r,t)  j  +  C^7x7xu®(f,t)  =  4nq(r,t)  (l-5) 


where  4TTq  « 

Po 

This  is  the  standard  form  for  the  inhomogeneous  vector  wave  equation.ii/  From 
now  on  we  shall  be  interested  only  in  the  scattered  wave  sp  we  shall  drop  the 
superscript  s  until  further  notice.  We  introduce  the  time  Fourier  transform 
pairs 


=  J  ei<"^u(?,t)dt 


'  k  / 


e  u(r,a)Jau) 


€ 


e^‘®^q(r,  t  )dt 


q(r,t)  =  ^  f  e~^*^q(r,w 


)diu 


(1-6) 


Prom  Eq.  (l-5)  the  equation  for  the  Fourier  transform  is 


-  U) 


I)^u{r,«))  -  C^7  |7*u{r,(u)|  +  C^73C7Xu(r,u))  =  4TTq(r,uj) 


From  Eq.  (l-4)  and  Bq.  (l-6)  one  finds  that 


q[r,u))  = - <  e  6tu)+(UT^+uta 

4iPo  L 


4ip 


} 


f  -iCv^2).?  ^  ^  Kvkg).? 

—  <e  6 )  -  e  Stuu-ouj^+iug ^ 


(I-V) 


(1-8) 
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where 


6(a)  =  0  for  a  ^  0 
6(a)  =  o>  for  a  “  0 


eind 


6(a)da  = 


X  . 


Now  Eq..  (1-8)  is  not  quite  right.  We  are  interested  in  the  case 
where  the  primary  beams  are  well  collimated.  The  interaction  term  ^  ,  which 
is  a  product  of  the  anqplitudes  of  the  primary  beams,  will  be  zero  unless  we 
are  in  the  region  where  the  primary  beams  intersect.  We  shall  denote  both 
the  region  of  intersection  and  its  volume  by  V  .  Ihen  Eq.  (l-8)  is  valid 
only  Inside  V  and  q  is  zero  outside  V  . 


If  one  assiiines  that  fl(r,u))  decreases  at  least  as  fast  as  1  for 

large  r  ,  then  it  can  be  showniV  that  the  solution  for  Eq.  (l-7)  in  the 
infinite  region  can  be  written  as 


Tl(F,u)) 


G(r.r'a)5(f '  ,(i))dV 


(1-9) 


where  G  is  a  tensor  or  dyadic  operator, 
be  writtenii/ 


For  the  infinite  region,  G  can 


C? 


(I-IO) 


In  Eq.  (1-9)  the  term  arising  from  Cjj  will  be  longitudinal,  while 
the  term  arising  ftom  G^  will  be  transverse.  In  dyadic  notation  G/  and 
G^  are  given  byi!/  * 
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L  (1-12) 

where  «  =  ?-?•  (Figure  4)  and  R  =  1r|  and  I  is  the  unit  dyadic. 


We  now  choose  r  so  large  that 


a'/CjjR  »  1  >  lu/c^  R  »  1 


(1-13) 


for  all  r’  in  V  .  If  in  addition  j?'!  «  |r|  we  have 


_  A  _i  , 

R  ^  r«r*r 


(1-14) 


where  r  is  a  unit  vector  in  the  direction  of  r  ,  i.e.,  r  -  f/r  .  With 
these  assumptions  we  can  simplify  Eq.  (l-ll)  and  Eq*  (l-12).  Firstly,  we  keep 
only  the  leading  terms  in  (u)/Cjg)R  and  ((a/c^)R  .  Secondly,  we  replace  R 
by  r-f-f'  in  the  exponentials  and  R  by  r  elsewhere.  We  obtain 


G£(f,f '  ,w/Cjj) 


A 


rr 


Gt(f,f’,u>/Ct)  = 


!iL 

c+ 


A  . 

r*r  ^ 


r 


(1-15) 


(I-16J 
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Figure  4  -  Vector  Arx^ispaenx. 
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Mow  using  Eqs.  (1-15 j  and  (l-16)  in  Eq.  (l-lO)  we  can  find  G  and 
using  G  in  Eq.  (l-9)  with  q  given  by  Eq.  (l-8)  we  can  find  .  Then 

from  Eq.  (1-6)  we  can  get  u(f,t)  .  The  result  of  this  straightforward  but 
tedious  procedure  is 


^  (r-?)  f 

Ancfp^  ^  ' 

H 

j(  c,  (c,  -‘) 

.  f 

4nC|pQr 

sin  < 

,  !--(?. I- )f  r 

4nC^Pj^r  dy 

sin  i 

(1-17) 


The  first  and  second  terms  in  Eq.  (l-17)  are  longitudinal  waves  with  the 
5\jnimed  frequency  i;j3^+'jU2  and  the  difference  frequency  u)]^-a)2  ,  respectively. 
The  third  and  fourth  terms  are  transverse  waves  vdth  the  summed  and  difference 

frequencies ,  respectively, 

i 

do  not  vary  during  the  f '  integration.  Now  look  at  the  first  term  in 

Eq.  (1-17).  As  we  integrate  over  r'  the  integrand  oscillates  with  frequency 

determined  by  the  coefficient  of  r’  ,  i.e.,  ^ — £  r-k^-kg  in  this  case. 

In  general,  the  result  of  this  integration  will  depend  on  just  how  the  waves 
fit  into  the  region  V  .  As  we  increase  V  the  value  of  the  integral  will 
oscillate  between  fixed  limits,  unless  we  can  find  a  direction  i-g  for  which 


The  terms  in  t,nc  .arguments  of  the  sine  function  involving 


(§ 


.ki-k2 


=  0 


(1-18) 
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If  we  can  find  such  an  Fg  ,  the  integrand  becomes  constant  for 
?  =  fg  and  the  amplitude  of  the  scattered  rave  in  this  direction  becomes 
proportional  to  the  voliaoe  of  interaction  V  .  By  Increasing  the  volume  V 
one  can  indefinitely  increase  the  amplitude  of  the  scattered  rave  in  the  di¬ 
rection  .  In  any  other  direction  the  amplitude  does  not  increase  in¬ 
definitely  with  volume  but  oscillates.  For  the  proper  choice  of  experimental 
parameters  this  will,  lead  to  scattering  which  is  sharply  peaked  in  one  direc¬ 
tion  and  whose  aiiq;>litude  is  proportional  to  the  volume  of  Interaction.  We 
shall  call  this  part  of  the  wave  the  scattered  rave  as  opposed  to  the  rest  of 
the  outgoing  rave  which  has  the  character  of  a  diffracted  rave.  Equation 
(l-18)  will  be  called  the  resonance  condition.  We  have  so  far  considered  only 
the  first  term  in  Eq.  (l-17);  however,  the  character  of  the  remaining  three 
terms  is  the  same  as  that  of  the  first.  The  resonance  conditions  for  the 
remaining  three  terms  are,  respectively, 


^  =  0 

=  0 


A 


(kj^-kg)  =  0  . 


III.  IHTERACTION  CASEF 


Three  cases  of  Interaction  between  two  Intersecting  waves  must  be 
considered,  namely, 

1.  Two  transverse  raves. 

2.  Two  longitudinal  raves. 

3.  One  transverse  and  one  longitudinsil  wave. 

We  shall  consider  the  case  of  two  transverse  waves  in  detail  by  way  of  exanqile. 
The  results  for  all  three  cases  are  given  later  in  Table  II. 
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For  the  case  of  two  transverse  waves  we  have 


Ao'ki  ~  60*^2  =  0 


Vfe  want  first  to  see  if  the  resonant  condition  Eq.  (I-IS)  can  be  satisfied. 
We  must  have 


or 


=  +  kl  +  S^j^'kg 


+  ^“*1^ 

cf  cf  c| 


where 


is  the  angle  be^ween 


and  k2  -  The  above  leads  to 


cos 


In  order  that  this  equation  be  satisfied  we  must  have  the  right  side  less  than 
one  and  greater  than  minus  one.  This  leads  to  the  condition  that 


1-C+/Cj  1+C+/Cjj 


For  any  uj]^/uj2  this  range  we  can  choose  the  angle  9  between  the  primary 
wave  vectors  and  ^2  that  we  get  a  scattered  wave  (appearing  in  the 

direction  of  tc]^+£2)  •  examines  the  resonance  conditions  for  the  last 

three  terms  in  Eq,  (l*17)  it  turns  out  that  none  of  them  can  be  satisfied  for 
this  particular  choice  of  primary  waves.  For  the  scattered  wave  then^  ve  have 
from  Eq.  (l-17)  only  the  longitudinal  wave 
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(1-19) 


(l+.K)  K 

- 5 - ^  ^  • 

4tiCjjPo 


Using  the  fact  that  Aq^Ie^  *  BQ»k2  =  0  in  this  case  we  have  from 

Eq.  (l-4) 

r  =  -  i  (h+aA)  |(Xo'4)(‘^l+kl4)  +  (4-%)(l|+2ki.^)Ao 

+  (Ao*k2>(^?^2ki-k2)\J  -  i  (KW3+A/4+B)(A^.B^)(k^.k2)(k2+k^) 

-  i  (A/4+Bl{Ao-k2)(to-ki)(V^l)  . 


The  amplitude  of  the  scattered  wave  depends  on  the  polarization  of 
the^primary  waves.  If  one  primary  wave  is  polarized  peiTpendicular  to  the 

plane  and  the  other  is  polarized  in  the  plane,  then  I =  0 
and  there  is  no  scattered  wave  even  though  the  resonan’  undition  is  met.  If 
both  primary  waves  are  polarized  perpendicular  to  the  ^,^2  plane  the  scat- 
tered  wave  amplitude  is 


-W4)  [(a=41)(c^-l).a(a.l)(c2*l)] 


(K+U/3+A/4+B)  [c^'^CSc^-l)  (a+a2)+c2(c2-i)  (a^+l)] 


(1-20) 


where  a  is  the  frequency  ratio  ujg/ujji^  and  c  is  the  velocity  ratio  C^/C,  . 
If  both  primary  waves  are  polarized  in  the  ^1,^2  Plane  the  scattered  wave 
amplitude  is  similar  in  form  to  Eq.  (1-20).  In  fact,  the  scattered  wave  ampli¬ 
tudes  for  all  three  general  interaction  cases  are  of  similar  form.  Since  these 
amplitudes  may  be  obtained  from  straightforward  expansion  of  the  terms  in  Eq. 
(l“17)»  they  will  not  be  included  in  the  present  paper.  The  general  character 
of  Eq.  (l-20)  will  be  discussed  in  a  later  section. 

1^  to  consider  the  case  where  both  primary  waves  are  longi¬ 

tudinal,  and  the  case  where  one  is  longitudinal  and  one  transverse.  Conplete 
results  are  given  in  Ihble  H.  For  the  cose  of  two  longitudinal  waves,  the 
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IMTERACTIOM  CASEb  WHICH  PRODUCE  A  SCATTERED  WAVE 
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resonant  condition  can  be  satisfied  for  the  fourth  term  of  Eq.  (l-17)  only. 

For  the  remaining  case  the  resonant  conditions  can  be  satisfied  for  the  first, 
second,  and  fourth  terms.  In  this  latter  case  of  interaction  between  a  longi¬ 
tudinal  and  a  transverse  primary  wave,  the  amplitude  of  the  longitudinal  scat¬ 
tered  waves  vanish  when  Sq  Is  polarized  perpendicular  to  the  plane  of 
ki,k2  .  Hie  trsuisverse  scattered  wave  has  finite  amplitude  for  any  polariza¬ 
tion  of  Bq  . 

The  appearsuice,  in  Table  II,  of  scattered  waves  with  the  difference 
frequency  requires  some  cominent.  Treatment  of  the  interaction  problem 

has  been  on  a  strictly  classical  basis  in  this  paper.  The  classical  treatment 
is  an  approximation  to  the  correct  Quantum  Mechanical  treatment.  A  macroscopic 
plane  elastic  wave  consists  of  the  presence  of  a  very  large  number  of  phonons 
of  a  particular  (long)  wavelength  in  the  crystal.  With  this  in  mind  one  can 
do  the  following  calculations; 

1.  Write  down  the  phonon  Hamiltonieuri  for  the  crystal  including  the 
first  term  giving  phonon-phonon  interaction. 

2.  Assume  that  at  t  *  0  ,  there  are  a  large  number  of  phonons  with 

wave  vectors  and  k2  and  energies  ,  -fttup  • 

3.  Assume  that  first  order  time  dependent  perturbation  theory  is 
valid  and  compute  the  state  of  the  system  at  some  time  t^  later  than  t  =  0  . 

If  one  does  this  calculation  the  following  points  appear.  Due  to 
the  structure  of  the  phonon-phonon  interaction  there  is  conservation  of  the 
phonon  wave  vector  (no  lAnklapp  processes  are  possible  if  the  initial  wave¬ 
lengths  are  very  long).  That  is,  a  phonon  of  wave  vector  and  one  of 

wave  vector  k2  can  produce,  in  interaction,  only  phonons  of  wave  vector 
^3  =  second  point  is  that  if  tj^  is  not  chosen  too  small,  the 

perturbation  theory  gives  only  energy  conserving  transitions,  that  is,  we 
must  also  have  -603  =  hc)2+ha)i  .  These  conditions  on  the  wave  vector  and 
energy  taken  together  sire  equivalent  to  the  resonance  condition  from  Eq. 

(I-I8).  A  troublesome  point  is,  however,  that  this  calculation  seems  to  indi¬ 
cate  that  a  scattered  wave  with  the  difference  of  the  primary  frequencies  can¬ 
not  arise,  in  contrast  to  the  classical  calculation  vdiich  permits  a  scattered 
wave  of  either  the  sum  or  difference  frequencies.  It  is  therefore  not  clear 
to  us  Aether  or  not  the  difference  frequency  waves  will  actually  be  produced 
in  an  e;q)eriinent.  Since  the  validity  of  the  assumptions  made  in  both  calcula¬ 
tions  is  somewhat  of  an  open  question,  it  is  difficult  to  make  a  more  definite 
statement  at  the  present  time. 
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IV.  BEAM  WIDTH  OF  SCATTERED  WAVE 


The  scattered  wave  given  oy  Eq..  (l-17)  appears  in  the  form  of  a 
conical  beam  with  vertex  at  the  interaction  zone  and  maximum  intensity  along 
the  direction  of  the  vector  fg  .  In  ejqieriments  aimed  at  detecting  the  scat¬ 
tered  wave  it  would  be  advantageous  to  minimize  the  spread  of  this  beam.  The 
parameters  which  determine  the  angular  width  of  the  beam  may  be  identified  as 
follows . 


Consider  the  first  term  in  the  right-hand  side  of  Eq,  (l-17).  The 
amplitude  of  the  integral  is  equal  to  the  volume  of  interaction  V  when 
f  =  f  .  We  seek  here  the  vector  for  which  the  amplitude  first  becomes 

zero  for  a  minimum)  as  f  moves  away  from  .  This  occurs  when  one  full 
cycle  of  the  wave  is  fitted  into  V  ,  or  in  other  words,  when  the  coefficient 
of  r'  is  approximately  equal  2tt/j6  ,  I  being  a  length  characterizing  the 
volume  of  interaction.  It  follows  that  the  angular  width  of  the  beam  is 
proportional  to  where  Xg  is  the  wavelength  of  the  scattered  wave. 

Thus  the  scattered  beam  will  be  narrow  when  Xg  «  i  .  For  the  case  we  are 
considering  Xs  is  inversely  proportional  to  03^+ wg  so  the  width  of  the 
scattered  beam  may  be  made  small  by  using  large  primary  frequencies  and/or 
a  large  interaction  volume. 

When  the  other  terms  in  Eq.  (l-17)  are  considered,  the  same  propor¬ 
tionality  to  Xs/j£  is  obtained.  However,  for  the  second  and  fourth  terms, 

Xg  is  inversely  proportional  to  u)]^-a2  which  suggests  that  in  general  a 
scattered  wave  with  the  difference  frequency  will  be  more  spread  out. 


V.  NUMERICAL  EXAMPLE 


We  use  a  numerical  example  to  Illustrate  the  previous  results. 
Choosing  polystyrene  for  the  elastic  medium  we  have^ 


X  = 


\1,  — 


K  = 


2.09  X  10^®  dynes/cm^ 
1.30  X  10^*^  dynes/cm^ 
3.01  X  10^®  dynes/cm^ 
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A*  *  -1,00  X  10^  dynes/cm^ 

B  =  -8.3  X  10^0  dynes/cm^ 

C  =  -1.06  X  10^  dynes/cm^ 

=  1.056  gm/cm^ 

Taking  again  the  case  of  tvo  transverse  primary  waves  we  find  that  resonance 
can  be  obtained  for 


0.338  <  a  <  2.955 


Since  the  two  primary  waves  are  the  same  type,  we  may  without  loss  of  gener 
ality  choose  loi  >  (U2  i*ich  restricts  a  to  the  range 


0.338  <  a  <  1 


The  smallest  angle  9  between  the  primary  waves  for  ^ich  a  scattered  (reso¬ 
nant)  wave  is  obtained  is  120.6  degrees  corresponding  to  a  *  1  .  As  a  de¬ 
creases,  9  approaches  180  degrees. 

If  both  primary  waves  are  poiarized  perpendicular  to  the  Tci  ,‘£o 
plane,  the  scattered  wave  amplitude  Eq,  (l-20)  is  largest  for  a  =  1  .  We 
then  obtain 


Amp.  (max.)  »  10.32  x  10’^  A^BoVcuf/r  cm. 


(when  Ao,  Bo,  and  r  are  in  centimeters,  V  in  cubic  centimeters  and  o, 
in  radlans/sec).  As  a  decreases,  the  amplitude  passes  through  a  minimum  at 
a  «  0.583  euid  we  find  that 


Amp. (min.)  =  7.45  x  lO-i®  AoBoV<aJ/r  cm. 


*  The  correlation  between  the  elastic  constants  i  ,  n  ,  and  n  given  in 
Ref.  3  and  A  ,  B  ,  and  C  used  here  is:  A  =  n  ,  B  .  m-^n  ,  C  .  i-m*in 
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Clearly,  the  amplitude  is  more  sensitive  to  the  primary  ’-ave  frequencies 
characterized  by  ,  than  to  the  frequency  ratio. 

Let  us  further  evaluate  the  interaction  in  polystyrene  of  two  10- 
megacycle  waves  with  a  volume  of  intersection  equal  to  1  cm5  If  the  displace¬ 
ment  amplitude  of  the  interacting  waves  is  approximately  10“10  cm.  and  we  let 
r  =  10  cm.,  then  the  displacement  amplitude  of  the  scattered  wave  Is  calculated 
to  be  approximately  10"^^ cm.  By  maintaining  the  same  volume,  frequency,  and 
observation  distance,  we  see  that  the  amplitude  of  the  scattered  wave  varies 
as  the  product  of  the  primary  wave  amplitudes.  Thus,  the  difference  between 
the  amplitudes  of  the  primary  waves  and  the  scattered  wave  decreases  as 
and  Bq  get  larger.  If  we  use  the  unreasonably  large  value  of  10"®  cm.  for 
Aq  and  Bq  (this  corresponds  to  a  strain  amplitude  of  about  10"^  in  the  above 
example),  then  we  find  that  the  amplitudes  of  the  intersecting  and  scattered 
waves  are  of  the  same  order  of  magnitude. 


VI.  CONCLUSIONS 


Two  intersecting  plane  elastic  waves  produce  a  scattered  wave  when 
the  resonant  condition  (Eq.  (l-18))  is  satisfied.  The  resonant  condition  is 
a  function  of  the  ratio  of  the  primary  frequencies  but  is  independent  of  the 
absolute  frequencies.  The  scattered  wave  appears  in  the  approximate  form  of 
a  conical  beam  emanating  from  the  volume  of  interaction  and  has  maximum  in¬ 
tensity  along  the  direction  fg  defined  by  the  resonant  condition.  The  width 
of  the  beam  is  proportional  to  Xg/i  \diere  Xg  is  the  wavelength  of  the 
scattered  wave  and  Z  is  a  length  characterizing  the  volume  of  interaction 
V  .  The  maximum  intensity  of  the  scattered  wave  is  greatest  for  high  primary 
frequencies,  large  primary  wave  amplitudes,  and  large  interaction  volume. 

In  an  experiment  aimed  at  detecting  the  scattered  wave,  the  distance 
L  from  the  interaction  zone  to  the  point  of  observation  should  be  large  com¬ 
pared  with  Z  .  Therefore,  optimxim  experimental  conditions  art 


L  Z  Xg 
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